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Abstract 


In this paper, we find radius of stability for all decision making units, with interval 
data. In this approach, organization classification remains unchanged under perturba- 
tions of the interval data. Some numerical examples for illustration are given. 
© 2003 Elsevier Inc. All rights reserved. 


Keywords: Data envelopment analysis; Interval data; Sensitivity and radius stability analysis 


1. Introduction 


Data envelopment analysis (DEA), introduced by Charnes et al. [1] (CCR) 
and extended by Banker et al. [2] (BCC), is a useful method to evaluate relative 
efficiency of multiple-inputs and multiple-outputs units based on observed 
data. Its goal is to classify the decision making unit (DMU) into two classes: 
efficient or inefficient ones. However, uncertainty such as a measurement error 
should be incorporated in observed data. This indicates the necessity to assess 
the sensitivity of classifications in DEA. In the recent years, many DEA 
researchers have studied the sensitivity of efficiency and inefficiency classifi- 
cations with respect to perturbations in data (see [3-6]). The original DEA 
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models assume that inputs and outputs are measured by exact values on a ratio 
scale. 

Recently, Cooper et al. [7] addressed the problem of imprecise data in DEA 
in its general form. Imprecise data means that some data are known only to the 
extent that the true values lie within prescribed bounds while other data are 
known only to satisfy certain ordinal relations. In this paper, we discuss a 
technique for assessing the sensitivity of efficiency and inefficiency classifica- 
tions in DEA with interval data. 

The focus is on the stability of DMUs that have interval data. Having 
identified efficient and inefficient DMUs in a DEA analysis with interval data, 
one may want to know how sensitive these identifications are to possible 
variations in data. 

In this paper, a modified CCR model is suggested to sensitive the DMUs 
with interval data. We develop several linear programming formulations for 
investigating radius of stability for all DMUs with interval data. The possible 
data perturbations for preserving the DMUs classifications are computed from 
the optimal values. 

Interval DEA models are extended to interval data [9]. Then, interval 
DEA for interval data can be extended to fuzzy data as well, since the level 
sets of fuzzy data are interval data. Therefore, fuzzy efficiency for fuzzy 
data can be obtained by interval DEA through the resolution identity 
[10]. 

The current article proceeds as follows: In Section 2, we review DEA 
models for dealing with interval data. Then, on the basis of these models, in 
Section 3, we propose some models for determining radius of stability 
for DMUs. In Section 4, the sensitivity and stability analysis methods to 
several data sets are introduced. A concluding section summarizes our main 
results. 


2. Preliminaries 


Suppose we have n DMUs, where each DMU,; (j € J, J = {1,...,}) con- 
sumes m inputs Xj! = (x1;,--.,Xmj) to produce s outputs Y' = (y1,---,¥). We 
assume Y,; > 0, X; > 0, ¥; 40, X;40. Unlike the original DEA model, we 
assume further that the levels of inputs and outputs are not known exactly, the 
true input and output data known to lie within bounded intervals, i.e. 
xi € [xy,x;;] and y,; € [yk,y5] with upper and lower bounds of the intervals 
given as constants and assumed strictly positive. In this case, the efficiency 
can be an interval. The upper limit of interval efficiency is obtained from 
the optimistic viewpoint and the lower limit is obtained from the pessimis- 
tic viewpoint. The following model provides such an upper bound for DMU,, 


[8]: 
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Ss 


=, U 
h ,, = max Se UY, 


r=1 


S.t. 3 vx, = 


i=1 

AY m ] 
Sowah, — Dooml, <0 se 
r=1 i=l 

Ss m 
So uy — d_ oxy <0, J=1,....0, JF jos 

r=1 i=l 


u,v Se, r=l,...,s, i=1,...,m. 
where ¢€ is a non-Archimedean infinitesimal and original variables w,,...,u; and 
Vj,---,Um are weights for outputs and inputs, respectively, which to be esti- 


msied We denote by hi the efficiency score attained by DMU,, in (1). 
The model below provides a lower bound of the efficiency score for DMU,,: 


AY 


a cE} 
h ;, = max S- UY, 


> U 
S.t. UX; = 1, 
i=1 


The efficiency hi. attained by DMU,, in (2) serves as a lower bound of its 
possible efficiency scores. Models (1) and (2) provide each DMU with a 
bounded interval (hi, sf ] in which its possible efficiency scores lie from the 
worst to the best case. Considering (1) and (2), it is evident that hi. < Ae: 

On the basis of the above efficiency score intervals, DMUs can be classified 


in three subsets as follows: 


fe opis : 
BE = {jet hi =}, 


Beste Jie SV cand. hy = 1}, 


BE ={j¢ J\hy <1}. 


466 G.R. Jahanshahloo et al. | Appl. Math. Comput. 156 (2004) 463-477 


3. Sensitivity analysis in DEA with interval data 


Suppose that DMUs are evaluated by (1) and (2) and are classified in E*T, 
E* and E-. Having identified efficient and inefficient DMUs in a DEA analysis, 
one may want to know how sensitive these identifications are to possible 
variations in the data. We determine “radius of stability” within which data 
variations will not alter a DMU’s classification from efficient to inefficient 
status or vice versa. We consider the radius of stability of DMU,, in three cases 
as follow. 


3.1. Radius of stability for DMU belonging to E** 


In this case, we assume that DMU,, is in E**, that is, hi, = 1. It is obvious 
that DMU,, remains in E** if its outputs increase or its inputs decrease. Our 
aim is to find the scalars « and f such that if we decrease upper bound of 
outputs and increase lower bound of inputs of DMU,, by « then h/’ = 1, also if 
we decrease lower bound of outputs and increase upper bound of inputs of 
DMU,, by f then he = = 1, ie. DMU,, remains in E**. 

We will calculate the upper- bounds of « and f, namely the ranges for « and 
f, such that the inputs and outputs of DMU,, are perturbed within these 
ranges and the efficiency of DMU,, is preserved in E**. It has been assumed 
that « and f are scalar and strictly positive. Here, we consider the following 
cases: 


(1) It is obvious if X/ — Ino, ¥7 + 1a then AY = 1 and if X/” — 1,8, Yi + 1,8 
then hie = =1; consequently, ‘DMU,, € E** (Note that 1 is a vector ith all 
components equal 1.). 

(2) If X¢ +10, YY —1a and X7 +1,,8, Yi —1,8, then it is possible for 
DMU,, not to ‘be i in the E*T, We are eonbenicd with finding the largest 
value fot a and f such that DMU,, € £*~. For this purpose, the following 
models are proposed: 


max «4 
s.t. Sou, + Sou =1, 
r=1 i=l 
S> Ur(Vr; a) 7 vilxi, + a) = 0, (3) 
r=1 i=l 
Sa eee 
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max fp 
r=l i=1 
Yo wr04, - B)~ Loaiel, +B) =0 (4 
Say! - ; vx, <0, fHl,....n, FA Jos 
r=1 i=l 
u,v; =e, r=l,...,s, i=1 m. 


Here, the constraint }~_, u, + 0", v; = 1 is added for normalization pur- 
pose. In the first glance, it seems that (3) and (4) are nonlinear, but the fol- 
lowing manipulation transforms the models into linear. 


Ss m 


S- u(y, =o) S- vi(x;;, +a) =0, 


r=1 i=1 


Ss m Ss m 

U L _ 
doh, — Dw, —%{ Dur + Do | =0, 
r=1 i=1 r=1 i=1 
Yunf, - Yet = 


r=1 


Similarly, we have 


By the above substitutions, the following linear models are obtained: 


max 


s.t. Se ei, 
r=1 i=l 


a Bet (5) 


r=1 i=1 


Souk — real <0, faleam Ade 


r=1 
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max fp 

s.t. yout vy, = 1, 
Yooh Set = 6) 
ye UY — Yoval fHlpseoryn, TF Jos 
uy,v; 2e, r=l,...,s, t=1,...,m. 


Assuming that «* and f* are optimal value of objective functions. Then, 
DMU,, with input vector X* +4 1,0* and output vector xe U _ 1.o* is in Ett, 


similarly DMU,, with input vector Xi’ +1,,6" and output vector Yi — 1B" isin 
Et. 


Theorem 1. For any a, where «> a", if xij, € [xi + Ina, x 


| @=1,...,m) and 
Yio © 5,1 ¥y, — Use] (7 = 1,.-.,5) then AY <1. 


Proof. In contrary, suppose there exists « such that « > a* and AY = 1. Since 


hi’ = 1, then there exists v;,u, (i= 1,...,m,r=1,...,s) such that 
dew t w= 1, 
r=1 i=1 


This means that « and v;,u, (i= 1,...,m, r= 1,...,s) are feasible solution to 
(3) and by optimality condition we have « < «*, which is a contradiction. 


Theorem 2. For any «, where « € [0,"], if xij, € [xjj, + Ina, x7] @=1,..-,m) 
and yrj, © [¥e, Yyj, — 1s0] (r= 1,...,8) then DMU;, € E* ie. he = =1. 


Proof. Suppose that there exists « € [0,«*] such that A < 1. Since hi’ < 1 then 
one can find 2 > 0 such that 


— et AP S (XG + Um) SOG, + Ano’) 
Seale J OA Ee She of ae 
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and strict inequality holds at least for one index. Considering model (3), this 
violates the assumption that DMU,, efficient one. 


Theorem 3. For any B, where B > B*, if xi, € [xt .x4 +18] G=1,...,m) and 


in i ijo?ifo 
Vrio © Wij, ~ 1,8, 5, (r= 1,...,8) then hy, <1. 


Proof. Considering the similarity between this theorem and Theorem | the 
proof is evident. 


By considering these facts the following theorem holds true as well. 


Theorem 4. For any B, where B € (0, B"], if xi, € [xt x2 +18] G@ = 1,...,m) 


ijo? vijo 


and y,;, € lye, —1sB, y9,] = 1,..-,5) then DMU,, € E** ie. hi, = 1. 
3.2. Radius of stability for DMU belonging to E* 


In this case, we assume that DMU,, is in E*, that is, hi. < land hy = 1. Our 
aim is to find the scalars « and f such that if we decrease upper bound of 
outputs and increase lower bound of inputs of DMU,, by « then h/’ = 1, also if 
we increase lower bound of outputs and decrease upper bound of inputs of 
DMU,, by f then Ai <1, ic. DMU,, remains in E*. 

We will calculate the upper bounds of « and f, namely the ranges for « and f 
such that the inputs and outputs of DMU,, perturbations within these ranges, 
preserve the efficiency of DMU,, in E*. Here, we consider the following cases: 


(1) If X; — 1,0, Yj’ + 1,0 then hj’ = 1 and if X,”, Y/ are fixed, then hi < 1, 
consequently DMU,, € E*. 

(2) If x” +18, ‘i —1,f then hi, < land eG yy are fixed then hi = 1, con- 
sequently DMU,, € E*. 

(3) We are concerned in finding the largest values of « and f such that 
DMU,, € E*. 


For this purpose the following models will be considered: 


max 0 
S.t. Sou + 54 =1, 
r=1 i=l 
dw, — 2) — Do oleh, +0) =9, (7) 
r=] i=1 
So uy — d_ vxt <0, J=1, fn, TF Jos 
r=1 i=l 
Up 0p 26. #= Vico, 8) TH 1.4m: 


470 


sup 
s.t. yey oe 1, 
r=1 i=1 
S- ur(Vi, + B)- ys vi(x;,, — 
r=1 i=1 
Dt 2 So 
u,,vu; =e, r=l,...,s, i=1 
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It seems that (7) and (8) are nonlinear, with similar manipulation the 


following linear models are obtained. 


max a 
Ss m 
S.t. ) uy; + ) vj = 1, 
r=1 i=1 
Ss m 
U a 
Yok, — road, = 
r=1 i=1 


r=1 i=l 
u,,v; >, r=l,...,s, i=1 
sup f 
AY m 
s.t be uy, + > v; = 1, 
r=1 i=1 
AY m 
U 
De Min — vy, +B <0, 
r=1 i=1 
S m 
UY; — SJ ox, <0, J=\1, 
r=1 i=l 


Instead of solving (10), the following problem is solved. For optimal value 
of , which is called £;, we have hi = 1, now we will prove that for all value of 


Rah hel: 
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min f, 


s.t. Sot 3c =, 
Sout, — Soeal, +B: = 0, (11) 


Assume «* and f} are optimal value of objective functions. Then DMU,, 
with input vector X/ +1, and output vector Y,/ —1,o* is in E*, similarly 
DMU,, with input vector x; —1,/ (6 < B}) and output vector ye +1,8 
(6 < B;) is in E*. The following theorems assert some concept concerning the 
fact stated above and some other result which are needed for our purpose. 


Theorem 5. For any a, where a>ot, if xy, € (xt +1nc,x¥] and 
Mis = WepIG, — Weel: then hy <1, 


Proof. Proof is evident. 


Theorem 6. For any «, where «€[0,"], if xy, € [xj + Ine,x;,] and 
Vein © ema — 1,0], then DMU,, € E* i.e. hi = 1; 


Proof. Proof is evident. 


Theorem 7. For any f, where B<B, if xi, © [ti xj, —InB] and 
rio © VG, + 158,97], then DMU;, € Et, ie. he <1. 


Proof. Suppose there exists # such that f < f; and hi; = 1. Since hi = 1, then 


there exists v;, u. G=1,...,m,r=1,...,s) such that 
Sat at 
r=1 i=1 


Ss m 


S-urlyt, +B) — D<v(x¥, — B) =0, 


r=1 i=1 


De — DP <0, j=1,...,n, JF Jos 


U,, Vj z €, 
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This means that f and v,, u, (i = 1,...,m,r=1,...,s) are feasible solutions 
to (11) and by optimality condition we have f > f; which is a contradic- 
tion. 


3.3. Radius of stability of DMU belonging to E~ 


In this case, we assume that DMU,, is in E~. Our aim is to find the scalars « 
and f such that if we decrease lower bound of inputs and increase upper bound 
of outputs of DMU,, by « then hi <1, also if we decrease upper bound of 
inputs and increase lower bound of outputs of DMU,, by f then hi < 1, Le. 
DMU,, remains in E-. 

We will calculate the upper bounds of « and f, namely the ranges for « and f 
such that the inputs and outputs of DMU,, are perturbed within these ranges 
and its inefficiency score is not changed. Now we consider following cases: 


(1) If AG, +10, Ke —1,« then hi <1 and if Xe +18, Ye —1,f then and 
hi, <1, consequently DMU,, € E-. 

(2) We are concerned with finding the largest value for « and f such that 
DMU,, € E-. 


We consider the following models: 


sup o 
S.t. Sou +o =1, 
r=1 i=1 
dour, +2) — Dwi, — 2) <0, (12) 
r=1 i=l 
So uye — So vxt! <0, J=1, A, TF Jos 
r=1 i=1 
u,v; >e, r=1,...,s, i=1,...,m. 
sup f 
S.t. Sou, + Sou = 1, 
r=1 i=1 
So u(y, + B) — So v(x, — B) <0, (13) 
r=1 i=1 
So uyy — vj i; <9, j=1, n, JF Jo, 
r=1 i=l 
u,v; 2€, r=l,...,s, i=1,...,m. 
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Similarly, the above mentioned model can be transformed to following linear 
models: 


sup « 
s.t. Sy owe 
r=1 i=1 
wal, Dewh, +2 <0, (14) 
r=1 i=l 
Daly Yoal <0, PH lyevign, J #63 
Uu,,0; 2€, r=l,...,s, i=1,...,m. 
sup fp 
s.t. ey ek 


Since the second constraints in (14) and (15) are strict inequality, then these 
models are nonlinear programming and the following linear programming 
problems are suggested: 


min 
s.t. ioe 
r=1 i=l 
use a5 ee +a, =0, (16) 
r=1 i=1 
Sus — Sond <0, fel, nceys JF Iai 
r=1 i=l 


Uu,,v; 2€, r=l,... 
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min fp, 
S.t. Sas ae 5 
r= i=l 
DMTin — DMG, + Br = 9 (17) 
Soa DLs j=1,...,2, JF Jos 
uy,¥; Se, r=1,...,5, i= 1,...,m. 


Assume a; and f; are optimal value of objective functions. Then DMU,, 
with input vector X/ — 1,0 (« < a) and output vector Y7 + Io (% < a7) is in 
E-, similarly DMU,, with input vector xX — 1, (6 < f;) and output vector 
¥ +1,f (6 < B}) is in E-. 


Theorem 8. For any «a, where a<aj, if xi, € ea - 1,0, x% ] and 
Vein © (Wij. Yq, + Used, then DMUj, € E ie. hi) < 1, where a; is the optimal value 


to (16). 


Proof. Considering the similarity between this theorem and Theorem 7 the 
proof is omitted. 


Theorem 9. For any f, where B < B}, if xi, © [X%j,.%4, — Im] @ = 1,...,m) and 


rj. © WG, +198, ¥5] 7 = 1,-..,8) then DMU,, € EW ie. hi, < 1, where Bi is the 
optimal value to (17). 


Proof. Proof is evident. 


Table 1 
Data of numerical example 
DMU, Input Output Efficiency unit Classification 
xy ty yyy yyy hy 
1 12 15 0.21 0.48 138 144 21 22 0.224 1 Et 
2 10 17 O11 0.7 143. 159 28 35 0.227 1 Et 
3 4 12 0.16 0.35 157 198 21 29 0.823 1 Et 
4 19 22 0.12 0.19 158 181 21 25 0.445 0.907 E 
5 14 15 0.06 0.09 157 161 28 40 1 1 Ee 


Table 2 


The sensitivity analysis result for DMUs 


DMU, Input Output Efficiency 
My xy X%y xy My My Yay Yay hy hy 
5 19.63626 15.00872 5.69626 0.09872 156.99128 155.36374 27.99128 34.36374 1 1 
Table 3 
The values of «* and f;, for DMU; 
DMUs ot By 
1 0.12445 0.42843 
2 5.69766 0.64648 
3 13.61625 0.27068 
Table 4 
The sensitivity analysis results for DMU; 
DMU;, _ Input Output Efficiency 
xi xi Xj x3) a vy Vay vay hy hy 
1 12.12445 14.5719 0.33445 0.0519 138.4281 143.87555 21.4281 21.87555 0.994 1 
2 15.69766 16.355 5.79766 0.055 143.645 153.30234 28.645 29.30234 0.973 1 
3 17.61625 11.74 13.77625 0.09 157.26 184.38375 21.26 15.38375 0.950 1 
Table 5 
The sensitivity analysis result for DMU4 
DMU; Input Output Efficiency 
J #4, 4 g a i i i if 
4 18.9839 21.87 0.1039 0.06 158.13 181.0161 21.13 25.0161 0.982 0.999 


LLE-E9b (FOOT) 9ST “IndwoD yiny jddy | Jp 1a oojyoysuvyye “YD 


SLY 


476 G.R. Jahanshahloo et al. | Appl. Math. Comput. 156 (2004) 463-477 


4. Numerical example 


Consider the interval data setting of Table | contain 5 DMUs with 2 inputs 
and 2 outputs and the efficiency scores obtained by applying (1) and (2) [8]. 

Considering the data from Table 1 and using models in Section 3, we cal- 
culus the radius of stability for all DMUs as the following: 


(i) Radius of stability E** 

For determining radius of stability of DMUs, we apply (5) and (6) (because 
DMUs is in E**). In this case, «* = 5.63626 and f° = 0.00872. If x5, for i = 1,2 
are increased by «* and y%, for r = 1,2 are decreased by «*, and apply (1), then 
hY = 1 and DMUs remains in £**. If xf, for i= 1,2 are increased by f*, and 
y;, for r = 1,2 are decreased by f", and apply (2), then Af = 1, consequently 
DMU,; remains in Ett. 

Table 2, reports the sensitivity analysis result for DMUs. 

(ii) Radius of stability E* 

By applying (9) and (11) (because DMU’s 1,2 and 3 are in E*) the data is as 
shown in Table 3. 

If xj, for i= 1,2 and j = 1,2,3 are increased by «; and yi, for r= 1,2 and 
j = 1,2,3 are decreased by «;, and apply (1), then hi = 1 and DMU, remains 
in Et, Lites for i = 1,2 and j = 1,2,3 are decreased by B; (B; < B;) and Vin for 
r= 1,2 and j= 1,2,3 are increased by f; (8; < B;), and apply (2), then hi < 1, 
consequently DMU; remains in E*. 

Table 4 reports the sensitivity analysis result for the DMU; where j = 1, 2, 3. 
In this case we assume f, = 0.4281, 6, = 0.645 and f; = 0.26. 

(iti) Radius of stability E~ 

For determining radius of stability of DMU4, we apply (16) and (17) (be- 
cause DMUs is in E~). In this case, «} = 0.01624 and f; = 0.1311. If x4, for 
i= 1,2 are decreased by « (a < af) and y¥, for r=1,2 are increased by « 
(a < «;) and apply (1) then AY < 1 and DMU, remains in E~. If x4, for i = 1,2 
are decreased by f (B < f;) and y4,, for r = 1,2 are increased by f (f < B;), and 
apply (2) then hf < 1 consequently DMU, remains in E~. 


Table 5 reports the sensitivity analysis result for the DMU4. In this case, we 
assume o = 0.0161 and 6 = 0.13. 


5. Conclusion 


In this paper, we proposed a method for finding the largest region for 
preserving classification of DMU with interval data. For this purpose, the 
classical CCR model is modified. Finding largest region which preserves the 
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classification, for each DMU a radius of stability is found. Such a work with 
fuzzy data may be done and the model may be extended. 
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